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Abstract 

We discuss the Hamiltonian for a nonrelativistic electron with spin in the 
presence of an abehan magnetic monopole and note that it is not self-adjoint 
in the lowest two angular momentum modes. We then use von Neumann's 
theory of self-adjoint extensions to construct a self-adjoint operator with the 
same functional form. In general, this operator will have eigenstates in which 
the lowest two angular momentum modes mix, thereby removing conservation 
of angular momentum. However, consistency with the solutions of the Dirac 
equation limits the possibilities such that conservation of angular momentum 
is restored. Because the same effect occurs for a spinless particle with a 
sufficiently attractive inverse square potential, we also study this system. We 
use this simpler Hamiltonian to compare the eigenfunctions corresponding to a 
particular self-adjoint extension with the eigenfunctions satisfying a boundary 
condition consistent with probability conservation. 

I. INTRODUCTION 

In this article, we first examine the Pauli Equation for an electron in the field of a 
magnetic monopole. This equation has appeared in the literature before, and it is well known 
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that an extension is needed to make the Hamiltonian self-adjoint in the j = sector [|l|. 
What seems to have gone unnoticed is that, for eg = ^, the domain to be extended includes 
the J = 1 sector as well. With the inclusion of this sector, the structure of the extensions 
becomes richer, and the number of parameters required to describe them jumps from 1 to 16. 
While the parameters may be chosen to be consistent with angular momentum conservation, 
this is not required: a pure incoming s-wave can come out with a p-wave component, even 
though the functional form of the Hamiltonian is spherically symmetric. However, if we 
require our states to match the states of the Dirac equation in the nonrelativistic limit, 
we will only have 1 free parameter, and angular momentum will be conserved. To better 
understand the effect of the extension parameters and their relation to angular momentum 
conservation, we also consider a simpler Hamiltonian of the form 

H = -—V^--^, (1) 

where c is an arbitrary constant. (Spin is not essential to this discussion and is omitted.) 
This Hamiltonian has all the essential features of the monopole Hamiltonian, so we can use 
it to investigate how the extension parameters arise. To accomplish this, we look at this 
Hamiltonian in 3-space minus a sphere of radius vq around the origin. We compare the 
results of imposing a boundary condition consistent with probability conservation with the 
results of creating a self-adjoint extension. Finally, we compare the case of a nonzero radius 
tq with that of a zero radius tq. 

II. PAULI EQUATION 

Working in units where the speed of light and Planck's constant are both equal to one, 
the Hamiltonian for an electron (with spin) in an electromagnetic field is 

H = — -—a-B, (2) 
2/^ 2/x ^ ' 

where vf is the kinematic momentum p — eA, A is the vector potential, and B is the magnetic 

field. For a point magnetic monopole of strength g sitting at the origin 0, we have 



/=^il^^0 (3) 
r sin 

(4) 

E- — - "^^-^ f5l 
^ " 2^ 2/ir2 ' 

where we have chosen a particular gauge to determine A. Using L = r ^ -k — egf, the 



Hamiltonian can be rewritten as 



Here L is the orbital angular momentum satisfying [U, V^] = ie^^^V^ for any spin- 
independent vector operator. Note that [U,H] ^ 0. The angular momentum operator 
that commutes with H is the total angular momentum J = L + ^a. In appendix I, is 
shown to simplify to 

1 1 92 K{K + 1) 

H = r H (7) 

2/i r dr'^ 2/ir^ 

with the further definition 

K=-l-rxn-a (8) 

p. We can find simultaneous eigenstates of K, J^, and Jz involving only angular variables. 
Such "monopole harmonics" will be denoted by fi^ml^', 0) and satisfy 

-^^Km I^^Km (9) 
J^n^n,=j{j + l)^Km (10) 
Jz^Km TIT'^Km (H) 

Jdn\n^^{e,<j))\^ = i (12) 

where j takes on the values 0, 1, 2, . . . and k, is related to j by 



K = ±^j{j + 1) (13) 
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for eg 



(More generally, k = ±J(j + ^y — e^g^ with j 
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See 



appendix I.) Note that for j > (or j > eg — }^\xy the general case), there are two sets 
of m-multiplets for each value of j, one corresponding to k > 0, and one corresponding to 
K < 0. In terms of Bessel functions and the monopole harmonics, the > solutions to 
the eigenvalue equation Hil) = Eip, are 



KmE 
N 



^KmE 

o 



(14) 
(15) 



and the E <Q solutions are 



^ B 



f 2 



(16) 



Vk 



(There is another set of < solutions which has ly^ instead of K^^, but they grow 
exponentially at large distances and need not be considered.) At this point we mention 
that any dependence on the gauge is contained entirely in the form of the ^nm- The radial 
part of these solutions and the eigenvalues of the operators are left invariant under a guage 
transformation. The set ([T^ ) (N for Nonsingular) of solutions vanishes at the origin, while 
the sets (^) (S for Singular) and (|16D (B for bound) are singular at the origin. Although 
the are not square integrable over all space, they can be (5-function normalized and are 
square integrable over any finite region. However, the \E's and b are only normalizable 
(5-function normalizable in the case of the s) when u,^ < 1 <^ k = 0, — -\/2 (i.e. for 
the j = singlet and one j = 1 triplet). If > 1, the and are not square 
integrable over any region containing the origin. (The cutoff z/^ = 1 is equivalent to a 
coefficient of k[k + 1) = | for the ^ term in (|^). This coefficient is a general cutoff for a 
^ potential |0| [||.) It is tempting to let the Hamiltonian operator (|^) act on all linear 
combinations of the normalizable (including (5-function normalizable) solutions. However, 
this is incompatible with the Hermiticity condition {Hcf)^ ip) = (0, Hip) for all 0, ip in the 



domain of H . Furthermore, we seek a Hamiltonian which is not only Hermitian, but also 
self-adjoint, for only then are its eigenfunctions complete. At this point, we need to become 
more precise in our usage of the term operator, from now on including the domain of an 
operator as part of the operator's definition. Precise definitions of Hermitian and self-adjoint 
will be employed, and can be found in appendix II. To start with, let Hi be the operator 
given functionally by in (^, and defined on all functions ip in the Hilbert space such that 
Hip is in the Hilbert space. We expect the Hamiltonian to be Hermitian; however, 

1 f.,r.(dr 



(Hiip, 0) = (^, H,<P) + lim - J r^dn ^-^<P - ^*^J (17) 

from integration by parts, so Hi is not Hermitian. Next, define H2 to be the operator 
identical to Hi except that the domain is further restricted to 

'dip* 



E dom(ifi 



lim J r^dn I - iP*-£ 1=0 for all G dom(i7i) I . (18) 



By comparing ( [T7| ) with (|T8|), if 2 is seen to be Hermitian. In other words, the domain of H2 
is the set of in the domain of Hi such that 

(V^,if20) = (iflV^,0) (19) 

for all ip in the domain of Hi, which means that if] = -^i- (-^1 is the adjoint of ii2.) Since 
the domain of ifi is larger than that of ii2, ii2 is not self-adjoint; however, the domain 
of ii2 can be extended through a method of von Neumann to create an operator that is 
self-adjoint. According to the von Neumann theory of self-adjoint extensions, we need to 
look at the number n+ of normalizable solutions to the equation iii0+ = and the 

number n_ of normalizable solutions to the equation Hi(p_ = —ifi(p_. (Note that the use of 
fi is arbitrary and chosen only to provide the correct units. Any positive real constant may 
be used instead.) We can index these solutions and denote them by (p"^^, where i ranges from 
1 to n±. If = n_ = n, then ii2 can be made self-adjoint by introducing the n vectors 
(p^ = (p\ + Uj(p-'_ into its domain, where U is an arbitrary unitary n x n matrix (with real 
parameters). Thus, a general element of the domain of the extended Hamiltonian, Hjj, is of 



the form Ci(j)'^ + (j), where (j) is in the domain of H2. (The reader can check that the extended 
Hamiltonian satisfies the self adjointness criterion of appendix II.) For the monopole, the 
domain of H2 consists of those functions in the domain of Hi that vanish at the origin at 
least as fast as r^/^. Then, the normalized solutions to Hi(j)± = ±in(f)± are 



V TT 

j = 0, K = 0, m = or j = 1, K = — v^, m = 0, ±1. 

There are 4 of each, so von Neumann's Theorem tells us that we need 4^ = 16 parameters 
to describe each extension. Given 16 parameters in the form of a unitary matrix Uj]^,, 
= 0, 1; m = — j, . . . , j; m' = —j', . . . , j'), the Hamiltonian can be made self-adjoint by 
introducing the 4 vectors 



^jm ^ ^jrn ^ ^^,^1"" (21) 



into its domain. (On a cautionary note: The superscripts jm on the (fP'"^ should be considered 
merely labels. The states (p^^ defined above and the '^■'^ defined below are not necessarily 
eigenstates of angular momentum.) The 0-^™ are not energy eigenstates; however, for each 
^-'"^ and for each positive energy eigenvalue E, there exists one energy eigenstate, ^l/^^, that 
differs from (p^"^ by an element in the domain of H2. When f/j/*^/ is diagonal, the '^■'^ are 
simultaneous eigenf unctions of J^, J2, and Hu, but in general this is not the case: The 
'^■'^ will be eigenfunctions of Hjj only, since the simultaneous eigenstates of i^i, J^, and 
corresponding to j = 0, 1 are not in the domain of Hjj. Similarly, not all of the eigenstates 
of and Jz will be eigenstates of Hu. In other words, a pure angular momentum eigenstate 
with j = 0, 1 will in time evolve as a superposition of states with mixed angular momenta. 
Thus, angular momentum is not conserved for general U . Now we can construct the energy 
eigenstates. Because the '^■'^ differ by 0-"" by an element in the domain of 7^2, which 
vanishes at the origin at least as fast as r^, we only need to consider the behavior of the 
solutions (|T^), (P^), and (|TB|) at the origin. To get a particular energy eigenstate, pick a 



value of energy E and a particular 0-^™. Then, look at its behavior and the behavior of the 
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above solutions for small r. Any linear combination of solutions whose small r behavior 
matches the small r behavior of the particular 0-''" (up to a part that vanishes at least as 
fast as r5) is an eigenvalue of the self-adjoint operator Hjj. For positive values of E, there 
is precisely one energy eigenstate for each 0-^'". However, most negative values of E fail 
to yield an eigenstate. When we try to match the solution to 0-''", we derive the following 



relation between the energy E and the diagonal matrix element U^^ : 



E = -fi 



1 + ^'^-^ 



1 



(22) 



where, as before, = 
K = for j = and k 



K + 



and K = ±yj{j + 1). (Specifically, we are only interested in 
for j = 1.) For the bound state to exist, we require that 



1. Since U is unitary, this 



E be real and negative. The reality condition requires 
requires the the row and column corresponding to j and m consist entirely of zeros, except 
for the diagonal element, which must be of the form e*^, where 9 is real. (Note that this also 
implies that the bound state must be an angular momentum eigenstate.) This allows us to 
simplify ( ^2]) to 

cos{TiVt^/2) + cos{9) 



E = -fi 



1 



1 + cos{e - TTiy^/2) 

So we see that we also have the additional condition on the diagonal element: 

cos{9) > —cos{T{Uf^/2). 



(23) 



(24) 



There exists one bound state for each such row; therefore, there can be anywhere from 
zero to four bound states, depending on the particular self-adjoint extension chosen. Now 
we compare our results with a similar treatment for the Dirac equation that has already 



appeared in the literature H |10|. We work with the Dirac Hamiltonian 



H = a ■ n + Pfj, 



(25) 



and in a basis where 
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cr 



a — 



a 



^ 0^ 







(26) 



and (J are the Pauli matrices. The advantage of this basis is that the Dirac spinor can 
be separated into upper and lower bispinors, where the lower bispinor is dropped in the 
nonrelativistic limit. When appropriately separated, the eigenvalue equation H%1) — Eip 
becomes 
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(27) 



U - E)^"'^ - i(dr + -)g''^^ 

r 

i{dr + i±^)/«"^^ + + E)^--^ 




0. 



(28) 
(29) 



Solving for y^™^^ we get the familiar solutions for E > (j, 



and ioT E < n 



fKmE 

Jn 



r ^Ju^{\r) 



X = jE^-fi^ ly. 



1 

^^2 



(30) 
(31) 



fKmE 

Jb 



X^yJiJ?- E\ v^^ 



1 



(32) 



These are the same solutions as for the Pauli equation, except that A has the different 
expression above. However, if we let £^ = + E' and identify E' with the nonrelativistic 



energy, we get A = 2\xE' + E' 



1/2 



^ \l\iE'\^^'^ in the nonrelativistic E' « /i limit, 
and we recover the solutions to the Pauli equation in this limit. One additional difference 
between the solutions to the Pauli and Dirac equations is that we now have an additional 
function g'^'^^ = —ii^dr + (1 + K.)r^^)f'^"^^/ {iJ, + E), which we require to be square integrable 
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over a finite region containing the origin. One may easily check for k = — -\/2 that g^^^ 
is not square integrable over this region, even though fg"^^ is. Thus, we no longer have 
the singular solution for k = —\pi as we did in the Pauli case above. When we look for 
self-adjoint extensions, we find that we no longer need an extension for the j = 1 sector; 
only the j = sector requires an extension. Thus, angular momentum modes may not mix; 
furthermore, we only need one parameter to specify the extension. From another point of 
view, consistency with the nonrelativistic limit of the Dirac equation requires us to fix 15 of 
the 16 parameters of the Pauli equation. 

III. GENERAL ^ POTENTIAL 

The Pauli equation in the presence of a magnetic monopole is just one example of a 
Hamiltonian where the choice of a self-adjoint extension can lead to non-conservation of 
angular momentum. Looking at the form of (0), we can see that the same essential behavior 
can be obtained from a spinless particle with a sufficiently attractive inverse square potential, 
as in Analysis of this new Hamiltonian is simple and analagous to the analysis of the 
magnetic monopole. Similarly, the arguments in this section can be modified to include 
the monopole Hamiltonian or any similar Hamiltonian that needs an extension. With this 
simpler Hamiltonian, we can investigate how the extension parameters arise. To do this, 
consider a space with a sphere of radius ro and centered around the origin removed. We 
impose conservation of probability at the boundary and seek a relation between the extension 
parameters and the boundary conditions. Now, consider a spinless particle governed by the 
Hamiltonian (|l]), whose functional form is given by (using L = —if x V) 



This Hamiltonian has appeared in the literature before. The solutions to Hip = Eip are 

^ivaE ^ r^\j^^^Xr)Y^^{p^) (34) 
mf'^ = r-W^X^r)Y^"'{Sl) (35) 



B 



r-2ir^,(Ar)r""(^]) 
\/7+i(j + l)-c, 



(36) 



where and \1>^^ are the solutions for > and \1>^^ is the only solution for < 0. 

'^^^^ and \E'^^ are still singular at the origin; nevertheless, because the sphere r < tq is no 
longer a part of our space, the singularity of solutions at r = is no longer important: All 
of the E <{] solutions can be normalized, and all of the E > {] solutions can be 5-function 
normalized. Now we wish to select those solutions that are consistent with probability 
conservation. Probability is conserved at the boundary r = ro if 



1*9, , ^ ,* 
or or 



0. 



(37) 



ro 



To ensure this, we impose the most general boundary condition at r = ro that is consistent 
with (|37|) by introducing a function (/^^(fi,^'): 



(Klgr^{Vt,Q!)i){rQ,Q!) 



(38) 



with the requirement that (^^.^(fi,^') is Hermitian, i.e.. 



(39) 



We allow the boundary condition to have a continuous dependence on rg through the explicit 
appearance of the subscript ro in (^) and (^). If gr^X^,VL') and ip{x) are expanded in 
spherical harmonics as 



grA^.n') = YUeA)C'^Y*^,[ 

^(f)=^'-(r)r,^(^,0). 



(40) 
(41) 



then ( |38| ) and (^) take on the matrix form 



dr 

I'm'lm* 



9 

9ro 



ro 

Iml'm' 



^''™'(ro) 



(42) 
(43) 
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Now, to construct the eigenfunctions of the Hamiltonian, we simply take the hnear combi- 
nations of the above solutions for a given energy E which are consistent with the boundary 
condition. The above procedure is self-contained and distinct from the von Neumann proce- 
dure It describes a boundary condition that restricts wavefunctions from the entire Hilbert 
space to a subspace on which probability is conserved at the origin. To make contact with 



the von Neumann procedure, we seek a relation between the boundary condition ( ^2|) and 
the unitary matrix that needs to be specified to apply the von Neumann theory. Instead of 
imposing the boundary condition (^) on all solutions to the eigenvalue equation, we can 
choose a particular self-adjoint extension Hu. To constuct Hu, start with Hi, the operator 
with the functional form of H in (|33D and domain consisting of functions ip in the Hilbert 
space such that Hip is in the Hilbert space. Then, create a new operator H2 as we did in (p!9D 
and look at the solutions to the equation Hi(f)± = ±ifi(j)±. This time, we obtain a solution 

= (P'^{r)Yim{9, 0) (no sum) (44) 

for each l,m. Since there are infinitely many /, m, and n_ are infinite, and we can create 
a self-adjoint operator Hu by extending the domain of H2 to include the infinite collection 
of vectors {0'™(x)}, where each vector is of the form 

0'-(f) = 0f (f) + UlZ'<P-"''{x) (45) 

and where U is an infinite- dimensional unitary matrix. (Again, the 0'"* are not angular 
momentum components in the sense of (^) and (0).) We now seek a relationship between 



f///™/ and 5''™''™''. Enforcing the hermiticity condition, {(j)'-^ , Huip) = {Hjjcj)'-"^ ^ip), for all -0 
in the domain of Hu, we have 

dm 



d'x<i)'"' {x)Huip{x) = / rf'x(iJ{70""(f))*0(f) 

r>ro Jr>ra 



(46) 



^ Im / rffi0-*(f) ^ = Im / rff] ^-^^iPix 
J or J or 

ro 

using the explicit functional form of the Hamiltonian, (|33|). Expanding ip as in (|^), we can 
take advantage of the orthogonality of the angular momentum harmonics to integrate them 
out, leaving only the radial part of the angular momentum components: 

11 



I'm' 



dr 



d 



ro 



I'm' 



(no sum on Im). 

(The "Im" of (^61) has been dropped since the phase of ijj is arbitrary.) We define 



ro 



a'-''™'(r) = (0f(r)4-, + f//C'0''™'(r))* (no sum). 
Note that a'™'''"'*(r) is the I'm' angular momentum component of ^'"^(a;), i.e., 

0'-(f)=a'™''™'*(r)llw(f^). 
Viewing a as a matrix in Im and I'm' which depends on r, we can write (^7|) as 

\ Iml'm' 



d 

dr 



i/;"^(r) 



a ^{r) — a(r) 



(47) 



(48) 



(49) 



(50) 



where the inverse and product inside the parentheses are a matrix inverse (but not a func- 
tional inverse with respect to r) and a matrix product. Since this is of the same form as 



ImVm' 
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Furthermore, if we take ip 
d 



Iml'm' 
ro 



a"m" 



(51) 



ro 



l"m"lm* 



in (pO|) (so that ip'-"^{r) = a (r)), we have 

Iml'm' 

J"m"l'm'* ■ 



d 



dr 



(r) = [a ^{r)—a{r) 



dr 



(r) 



d 



t 



d 



a {j')—a{r) = \a {r)—a{r) , 



(52) 
(53) 



dr J \ dr ^ 

which shows that g^^^'^' is Hermitian. Equation (|5TD is the desired link between the ex- 
tension and the boundary condition. Now that we have found it, we ask what happens in 
the limiting case tq = 0. To this effect, first note that, for a fixed value of c, most of the 
singular Bessel function solutions cease to be normalizable in this limit. The same holds for 
the 0f^(x): Since many of the 0f^(x) are no longer normalizable, any self-adjoint extension 
defined in this section which adds non- normalizable (f)^^{x) to the domain of Hi is no longer 
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valid. To determine which of the solutions and are still at least 

5-function normalizable, we require, as in the previous section, that the coefficient of the ^ 
term in the Hamiltonian be less than |. Then, our requirement becomes / < /^rit' where 

3 

W(W + l)-c=-. (54) 

Those extensions that remain valid have U diagonal for / > l^^j^f, with entries such that 
the linear combinations of and within the domain of are purely "^^^^ for 

/ > /^rit- ^^^^ point, the reader may object that we are including an infinite number 
of vectors in these extensions while the von Neumann indices are now finite. This is not 
a problem, since the with / > l^^j^f that we are including already exist within the 

domain of Returning to equation (^Tj), it follows that the only boundary conditions 
admissible in the ro — limit are those for which gl^^'"^' is diagonal, except possibly for the 
block with 1,1' < /grit as tq — 0. One should note that the relationship between Qrg and 
U becomes singular as tq — > due to the singularities of the 0^(r) in aj,^''™'. In general, 
finite entries in U lead to singular entries in Qro-^o- Thus, for tq = 0, it is more convenient to 
describe the domain of the Hamiltonian through U than through a boundary condition at 
the origin. However, if one asks for any physical description of the choice of extension, the 



formulation (42) is more valuable. It tells us, for instance how the radial flux leaving the 
boundary through the Im channel is related to the probability amplitudes at the boundary 
for each channel: 

^^-*(r)-^^'-(r) 
ar 



= E ^'""^ (ro)^'""''""'^''™' (ro) (no sum on Im) . (55) 

I'm' 

In summary, we can deal with the singularity at the origin by removing a small sphere of 
radius Tq from around the origin. When we do this, we must impose a boundary condition 
consistent with probability conservation, and we need a periodic function in two angular 
variables to describe this. With certain restrictions on the boundary condition when ro = 0, 
the angular momentum components of this function are in direct correspondence with the 
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elements of the unitary matrix required in the von Neumann theory. It is in this sense that 
the boundary condition is equivalent to a choice of self-adjoint extension, or alternatively, 
that the self-adjointness condition is equivalent to probability conservation at the boundary. 



APPENDIX I 



We start with a Hamiltonian of the form 



^ = (- ■ -)^(- ■ -) + • (56) 



Using —egr ■ a = (L — r x tt) ■ cr, we can rewrite the contents of the curly br aces as 

{} = L2 + L-<T-eV-(^X7r)-a (57) 
= (L + ]^df - 1 - (f X Tf) • a + i - e^g\ (58) 

Following the convention of |^ and 0, we define K = —1 — {r ^ -k) ■ a . For an eigenstate 
of .P (where J = L + this then becomes 

{} = J^ + ir + ^-eV- (59) 

The eigenvalues of will be of the form j{j + 1) for j = eg — ^,eg + ^, . . .. We will show 
below that the eigenvalues of K are 



«: = ±^(j + -)2-eV. (60) 

Given this relation between the eigenvalues of K and J^, the operator represented by the 
terms in curly braces has eigenvalues k{k,+ 1), where k, is the eigenvalue of K. This demon- 
strates the equivalence of the terms in (^ and (^. Next, the vector potential (0) for the 
magnetic monopole has no radial component, so n ■ f = p - r and f ■ vf = f-p, where p is just 
the mechanical momentum operator given by — iV. Thus, the first term of (P) is just the 
usual —^^^1^ term. To prove (|60|) , first write 

(r X vf) X (r X tt) = (L + egr) x (L + egr) = L x L -\ {L x r + fx L). (61) 
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Now, L satisfies [L*,V^-^] = ie^^^V^ for any vector operator tliat does not depend on spin 
coordinates, so tlie relations L x L = iL and Lxr + rxL = 2ir liold, and 

(r X 7?) X (r X tt) = i{L + 2egr). (62) 

If we now square tlie operator K, we see tliat 

= l + 2a ■ {fx n) + {a ■ {fx n)f (63) 
= 1 + 2(j ■ (f X tt) + (f X tt)^ + ia ■ [{f x f) x {fx tt)] (64) 

using the identity {a ■ A)^ = + ia ■ {A x A) . Making the substitution f x n = L + egr 
{{fx vf)^ = L? — e^g"^), this finally becomes 

K^ = l + a-L + L^- e^g"" = + ^ - e'g'- (65) 

Now, K commutes with since it is a total angular momentum scalar. Therefore, we can 
find eigenf unctions fi^m of K and with respective eigenvalues k and m. (Although we do 
not prove it here, these eigenfunctions are complete.) Then, from (|65|) , 

K^n^^ = ( J2 + 1 _ e^g^)n,^^ (66) 

which shows that fi^m is also an eigenfunction of J^. Replacing in (^) with its eigenvalue 
j{j + 1), we then arrive at 

«:^=j(j + l) + ^-eV = (j + ^)^-eV. (67) 

This proves (pDD, as long as we can show that k, takes on both positive and negative values. 
The latter follows from the anticommutation of K with {a ■ f): 

K{a -r) + {a ■f)K = 0, (68) 

as the reader may verify. Then, if KQ^m = i^^nm, 

K{a ■ r)n^rn = -{^ ■ r)Kn^rn = -^{a ■ r)n^rn (69) 
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(the m value is unchanged through multiphcation hj a ■ r). Furthermore, it can be shown 
that 

{a ■ r)Qoo = ^00- (70) 
So, from (^) and a phase convention that is consistent with ([70|) , 

fl-^m = (<? • r)Qi^m- (71) 



APPENDIX II 



The definitions of Hermitian and self-adjoint are [0: An operator H is Hermitian 
if its domain is dense, meaning every state in the Hilbert space can be arbitrarily well- 
approximated by states in the domain, and if 

(i/V^,0) = (^,i/0), (72) 

for every and in the domain of H. The adjoint of a densely defined operator H, denoted 
if^, is defined for any ip such that there is an t] for which 

{^,Hcl>) = ir],<P), (73) 

for all in the domain of H. In this event H^ip = 7]. H is self-adjoint ii H = H\ A crucial 
part of this definition is that the domains of H and W be equal. For a Hermitian if, the 
domain of W is always at least as big as the domain of if, and Wcf) = H(f) if is in the 
domain of H. 
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